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Exact Equilibria of a Stellar System in a Linearised Tidal Field 
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ABSTRACT 

We study the motion of stars in a star cluster which revolves in a circular orbit about its parent 
galaxy. The star cluster is modelled as an ellipsoid of uniform spatial density. We exhibit two 
2-parameter families of self-consistent equilibrium models in which the velocity at each point 
is confined to a line in velocity space. We exhibit the link between this problem and that of a 
uniform rotating ellipsoidal galaxy. With minimal adaptation, Freeman's bar models yield a 
third family. 
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1 INTRODUCTION 

In a star cluster moving in a circular orbit about the axis of sym- 
metry of a galaxy, a star is subject to the self-gravity of the cluster 
and the tidal field of the galaxy. If studied in a rotating, acceler- 
ating frame following the orbital motion of the cluster, it is also 
subject to centrifugal and Coriolis forces. If the spatial distribution 
of stars is uniform within an ellipsoid, and if the usual linear ap- 
proximation of the tidal field is adopted, the accelerations inside 
the ellipsoid are all linear in the spatial coordinates. There are then 
two norma l modes of oscillation in pl anes orthogonal to the axis 
of rotation. Fellhauer & Heggie ( 2005) studied this problem in the 
case when one of the normal frequencies is imaginary. Taking for 
convenience the case of an ellipsoid of revolution, they showed that 
it was possible to choose the axial ratio so that it was equal to that 
of the orbits in the remaining normal mode. By a very simple orbit 
superposition it was possible to construct a self-consistent model. 

This problem closely resembles that of constructing a model 
of a uniform, rotating ellipsoidal bar (Freeman 1966a,b,c; Hunter 
1974, 1975). This goal has been achieved in three cases: (i) ellip- 
tical cylinders, in which the semi-major axis parallel to the axis of 
rotation is infinite; (ii) so-called "balanced" systems, in which the 
centrifugal and gravitational forces are equal along the major axis; 
and (iii) disks with the surface density that is obtained by projecting 
a uniform ellipsoid onto one plane. 

In this paper we set out the connection between these two 
problems, and give a more systematic treatment of the models 
which Fellhauer & Heggie stumbled upon. In principle there are 
three dimensionless parameters: two axial ratios, and one parame- 
ter which measures the density of the system. The requirement of 
choosing one axial ratio appropriately reduces the family to a two- 
parameter family. Fellhauer & Heggie chose to work with axisym- 
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metric models, which simplifies consideration to a one-parameter 
family. In cases where both normal frequencies are real, however, 
it is clear that there is a possibility of a second two-parameter fam- 
ily. Finally, when this second eigenfrequency is zero, the problem 
is equivalent to that of Freeman's "balanced" bar. 

We begin in the next section with a summary of the equations 
of motion and their solution. Then section 3 surveys the models 
based on a single normal mode, and those which exploit the idea 
behind Freeman's bar model. The paper concludes with some fi- 
nal discussion, and an Appendix describes the models in terms of 
distribution functions. 



2 ORBITAL THEORY 

We use the rotating, accelerating frame in which the centre of the 
ellipsoidal star cluster is at the origin, the a>axis points to the centre 
of the galaxy, and the y-axis points in the direction of orbital motion 
of the system. (Fellhauer & Heggie took the y-axis in the opposite 
direction, but we switch in order to increase the resemblance with 
the theory as set out by Freeman.) Then the equations of motion are 

x + 2Qy + (k 2 - 4tt 2 )x = -A 2 x 



y - 29.x 

Z + V Z 



-By 



(1) 



where 9 is the angular velocity about the galaxy, k is the epicyclic 
frequency, v is the frequency of small motions orthogonal to the 
galactic orbit of the cluster (if the cluster itself is neglected), and 
A 2 ,B 2 ,C 2 are coefficients of x 2 /2, y 2 /2 and z 2 /2 in the expres- 
sion for the gravitational potential inside a uniform ellipsoid with 
semi-axes a, b, c. (They are functions of the density (p) and the ax- 
ial ^atios_fe/ ^c/cl) For a > b > c, these expressions are given 
in lFreemanl dl966bh . eqs.(4)-(12), which we would refer to hence- 
forth as FII.4-12, with corresponding abbreviations for equations 
in his tw o other papers. Equivalent f ormulae will be found in Ta- 
ble 2-1 of lBinnev & Tremaind fl987h . along with the special forms 
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corresponding to axisymmetrical and spherical cases. Thus 



x + 2Qy + Afx 
y - 2VLx + Bly 



z + Cfz 



where 

A\ - 

cl - 



k - AO. 2 + A 2 , 



B , and 

v 2 + C 2 



(2) 



(3) 
(4) 
(5) 



(cf. FI.8, FII.16). 

For a 7^ the solutions of these equations are 

x — A a sin(ai + e a ) + Ap sin(f3t + ep) 

y = k a A a cos(at + e a ) - kpAp cos(f3t + ep) 

z — A-y sml^yt + e 7 ) 



(6) 



(FI.9, FII.17c), where A a , Ap, A 7 are arbitrary constants, which 
we can take positive, £ a ,e,3,e 7 are arbitrary constants, 7 = Ci, 
a, f3 are the two roots for £ of the quartic 

C 4 - (A 2 + B'f + m 2 )f +A\Bl = 0, 



and 



k a — 
kp = 



A 2 2 

A\ — a 
2aQ, 

f3 2 -A\ 
2/3Q 



(7) 
(8) 

(9) 



(FI. 10-11). We adopt Freeman's convention that a < (5 when both 
frequencies are real; when one frequency is imaginary (a case not 
considered by Freeman) we denote the single real frequency by 
p. Then both k a and kp are positive. Note also that, for a given 
galactic potential, they are functions of p, b/a and c/a. 



3 THREE FAMILIES OF MODELS 
3.1 (3 models 

We first summarise the fa mily of self-consistent mo dels which 
generalises those found by iFellhauer & Heggiel d2005h . We take 
A a — in eqs.{6]l. The models are then built from the normal mode 
with frequency f3, and so we refer to them as /3-models. We can 
construct self-consistent models provided that kp — b/a. When 
the frequencies j3 and 7 are incommensurable, each orbit fills the 
curved surface of an elliptic cylinder. The upper end of this cylinder 
is the ellipse x 2 + y 2 /kp — A 2 p, z — A y , which lies on the surface 
oftheellipsoidxV^+j/V^+^/c 2 = 1 if A 2 /a 2 +A 2 /c 2 = 1. 
Thus in general 

4 ^Al max = c 2 (l-Al/a 2 ), (10) 

provided that A 2 p ^ a 2 . 

For given A/3,A-, the distribution of z is f(z\A-,) = 

= , 1 2 1 ^ A-,. Given in, we choose the distribution 
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(11) 



(For the case a — c, considered by Fellhauer & Heggie, this is 
equivalent to their eq.(19).) With this choice we see easily that z is 
uniformly distributed in \z j ^ A 1 , max ■ Finally, in order to generate 
a spatially uniform density, we choose the distribution of Ap so 
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Figure 1. Contours of equal density for the /3-models. The value against 
each contour line is GpjVi 1 



that f(Ap)dAp is proportional to the mass inside the cluster in the 
cylindrical volume between x 2 + y 2 /kp = A 2 , (Ap + dAp) 2 , i.e. 



(12) 



f(Ap) = ±Ap^l-Al/a\Ap < a 

(cf. Fellhauer & Heggie, eq.(18)). While these formulae are conve- 
nient for constructing a model, we also (section lA2t provide a more 
conventional definition in terms of a distribution function in phase 
space. 

Our purpose now is to sketch the range of models which sat- 
isfy the essential requirement kp — b/a. Since kp is a function of 
p, b/a and c/a, this equation may be solved for p as a function of 
b/a and c/a. This can be exhibited by plotting contours of p in the 
plane b/a, c/a. Results of a numerical survey are shown in Fig.l, 
where for definiteness we consider the case of an isothermal galaxy 
potential, in which k = \f2~Q,. 

Two limiting cases can be discerned: 

(i) Low density: If p = then f3 = k and kp = y/2, by eqs.Q, 
((9}. This is the limit at the extreme right of Fig.l. 

(ii) Axisymmetric models: If p is very large and a = b, we have 
A — B , ft ~ A + £1 and kp ~ 1, by the same pair of equations. 
This limit is reached at the extreme left of Fig. 1 . 

The models discussed by lFellhauer & Heggie lie on 

thelinec/a = 1 in Fig.l. (This subfamily is also depicted in Fig.4.) 
We now see that they are part of a two-parameter family. 



3.2 a models 

For sufficiently low densities, the value of a 2 is negative, and or- 
bital motions are unstable. For higher densities, however, a 2 ^ 0, 
and so we have another family of self-consistent models where now 
Ap = and k a = b/a. These we refer to as a models. Their dis- 
tribution in parameter space is shown in Fig. 2. Again there are ax- 
isymmetric models of arbitrarily high density. The construction of 
these models is exactly analogous to that of /3-models. Again the 
subfamily on which c/a — 1 is depicted in Fig.4. 



3.3 Freeman's Bars 

In F.II, Freeman described a series of three-dimensional ellipsoidal 
models of uniform density for the case in which A\ — 0. In the 
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Figure 3. Contours of equal density for Freeman's Bars, labelled as in Fig. 1 . 
The series ends at the near-vertical curve on the right, where kp = b/a. 



galactic context this was referred to as the condition for a "bal- 
anced" ellipsoid. For star clusters the corresponding interpretation 
is that of a system which is marginally stable against tidal disrup- 
tion, as a — 0. They appear to exist only for Gp/O 2 > 0.3 ap- 
proximately, and their distribution is shown in Fig. 3. 
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Figure 5. A typical orbit in the spherical Freeman model. 



For the case considered here, eqs.l|6^,b) are replaced by 

x — A a + A/3 cos((3t + ep) 
y = kpAg sin(/3t + e ) 



(13) 
(14) 



(FII.17, with some changes of notation), where now (3 2 — 4fi 2 + 
B 2 , k@ — 13/(20,), and the constant term A a represents the zero- 
frequency mode. The effect of this term is that epicycles can be 
placed with guiding centre at any point along the a>axis. This free- 
dom allows the construction of models with an axial ratio different 
from that of the epicycles. A little thought, however, shows that 
the axial ratio of the epicycles (kp) must exceed that of the cluster 
(b/a) in order that all parts of the cluster are accessible by epicy- 
cles lying entirely within it. This condition (which is best thought 
of as an upper limit to b/a) gives the near-vertical curve at about 
b — 1.2 in Fig[3j realisable models must lie to the left of this line. 
This curve represents the intersection of the families of Freeman 
and beta models, as can be seen in the cross section in FigU 

An attractive special case is the spherical model a — b = c. 
Then A = B = C = 4yrGp/3, and so eq.©, with Ai = 0, yields 



4ttG 



(AO 2 



. A typical orbit is shown in Fig[5] 



Figure 4. A cross section through all three families for c/a = 1. 



For Freeman models it was found convenient to generate ini- 
tial conditions from the distribution function, which is described in 
FIT. 



4 CONCLUSIONS AND FINAL REMARKS 

In this paper we have considered the problem of constructing an 
exact self-consistent equilibrium of a star cluster in a tidal field. 
The usual linear approximation of the tidal field is adopted, and we 
assume that the motion of the cluster in the galaxy is uniform and 
circular. 

We have found three distinct families of models, in all of 
which the cluster is an ellipsoid of uniform density. Each is a two- 
parameter family, the two parameters being the axial ratios of the 
ellipsoid. Additional free parameters are the angular frequency of 
the motion around the galaxy and the radius of the cluster. 
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All three models are built of epicyclic motions, of the kind 
which are familiar in galactic dynamics, though the details of 
the orbits are modified by the action of the cluster potential. In 
two families of models, on e of which was partially explored by 
iFellhauer & Heggiel d2005h . the space density is built up by epicy- 
cles whose guiding centre lies at the centre of the cluster and whose 
axial ratio (in the x, y directions) coincides with that of the ellip- 
soid. We call these a- and /3-models. In a third family the guiding 
centre of the epicycles may be displaced at various points along the 
x-axis, which points to the galactic centre. This family closely re- 
sembles Freeman 's family of "balanced" models for barred galaxies 
dFreemadl 196 6b). One of the families has models for all values of 
the space density (for suitable choice of axial ratio), while the other 
two exist only for densities above a certain minimum value. 

We have given enough details of the distribution functions as- 
sociated with these models to allow explicit iV-body realisations of 
them to be constructed. Code is available from the authors for the 
construction of the a- and /3-models, and for the spherical Freeman 
model. 

An interesting extension of the study of these models would 
be consideration of their dynamical evolution, under both collision- 
less and collisional processes. For a one-dimensional subfamily o f 
the /3-models this was considered by IFellhauer & Heggie I f2003) , 
with particular emphasis on low-density models. They found that 
the models are unstable, and that the time scale of the instability 
increases as the density decreases. Though they did not explore 
high-density models in such detail, their results indicated that the 
lifetime increases again at the highest end of the range they studied. 

The question arises whether it is possible to construct more 
general models along these lines. Though no further analytic mod- 
els have been found yet, it has been pointed out to us by J.P. Os- 
triker that this proble m seems well suited to Schwarzschild's orbit- 
superposition method jSchwarzschiidll 19791) . 
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APPENDIX A: DISTRIBUTION FUNCTIONS 
Al Integrals of the motion 

We invoke Jeans' Theorem and express the distribution function as 
a function of the isolating integrals of motion. As shown in FI and 
FII these are 

(\ 2 
■ . Pv 



f3k 2 _ &Pkp 



L: = -±A 9 
c, 12,2 1 .2 . 1 2„2 



where — = ak a + j3kp. These are in fact the energies of the three 
a _ 

normal modes, and the total energ}[jis E = E a + Ep + E 7 . 

By combining results in FI.41 and FJX27, we see that the con- 
dition that the envelope of this orbit just touches the surface of the 
model is J = 1, where, in general, 



J = 2o(k%-kl) 



E a 



+ 



ak a (a 2 k 2 - b 2 ) I3kfi{b 2 - a 2 k 2 ) J 7 



+ - 



•(A2) 



This expression is very useful for constructing distribution func- 
tions vanishing outside the surface of the model. When a is imag- 
inary only the /3-models exist, and then E a — 0. When a = we 
have the case of Freeman models, and then we find somewhat dif- 
ferent expressions for J and the quantities it contains (FII). In gen- 
eral, FI.12-19 are a useful resource of handy identities for checking 
some of the expressions in the following subsections. 

A2 a- and /3-models 

We consider here only the /3-models, as the distribution function for 
Q-models can then be constructed by obvious substitutions. When 
A a = E a — and b/a — k/3, eq. JA2b becomes simply 



J = 



2oE B 



+ ■ 



f3a 2 kp 7 2 c : 



(A3) 



Our search for the distribution function can be guided by 
eq.lll2l>. Substituting the first part of eq. JAlb ) into eq.d 10ft, we 

,2 2 ( 2(7 \ 

see that Ai ftmax = c II Ep I . From the first part of 

\ a pkp J 

eq. JAlb ). and eqs.QTJ and dA3t . it follows that the probability den- 



sity function of E^ is f{E 1 



-f(Aj) oc 



1 



I 2 A-/" " STZTj 



In fact 



it can be seen that this expression, with a ^-function to set E a = 
and with suitable normalisation, gives the full distribution function: 



f(x,y,z,x,y,z) 



pact 



S(E a 



2n 2r yc 



2aE ri 



2E-, 



(A4) 



f5a 2 kp 7 2 c 2 



It is easy to check that this yields the correct density, by integrating 
over the velocity: for the integration with respect to x and y we 
convert to suitable polar coordinates, which leads to an integral of 

f°° 2 

the form / 8(r )2rdr, and for this we adopt the value 1. 
Jo 

This paper has been typeset from a TjX/ KTpX file prepared by the 
author. 

1 in the rotating frame, so this is really the Jacobi integral 



